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1. Introduction 

Throughout this paper we denote hy B — {Bt,t^{)} a one-dimensional Brow- 
nian motion starting from defined on a probabihty space (fJ, !F , P) and by 
— {b[^^ = Bt + fit,t ^ 0} the corresponding Brownian motion with con- 
stant drift /i e R. In this and the second part of our survey (,51j)i we are 
concerned with the additive functional A^^^ — {aI''-*} defined by 



4"' = / exp(2B(^))ds, t ^ 0. 
Jo 



When /i = 0, we simply write At for aI'^K 

Such an exponential type functional plays an important role in several do- 
mains, e.g., mathematical finance, diffusion processes in random environments, 
stochastic analysis related to Brownian motions on hyperbolic spaces. This fact 
motivates detailed studies about this functional. 



*This is an original survey paper. 



312 



H. Matsumoto, M. Yor/ Exponential functionals of BM, I 313 

In this first part, we mainly consider the distribution of A[^^ at fixed time t, 
while some related results are also shown. It should be mentioned that, al- 
though A^'^^ is of a simple form, only integral representations for the density 
(of somewhat complicated forms!) are known by the results due to Alili-Gruet 
13, Comtet-Monthus-Yor [T3], Dufresne ,23^, Yor See also [3,111, |13|, EDI, 
ins and EI]. 

We mention the results obtained in dJ, [23 and |^ in some detail 
and we show some related identities, which are of independent interest. For 

example, the Bougerol identity sinh(i3t) ^ = ■* for fixed t and for another 
Brownian motion {Pu} independent of B and the Dufresne recursion relation, 
which relates the laws of A^'^^ and a['^^ (different drifts), will be mentioned. 
We also show some results about the moments of a[^^ and their asymptotics as 
t oo. 

The organization of this paper is as follows. 

In Section |5] we introduce the Hartman- Watson distribution after recalling 
some facts on Bessel processes. 

In Section 131 we show the Bougerol identity. 

In Section 0] we present several expressions for the density of Aj^-* at a fixed 
time t and at an independent exponential time. 

In Section we present the computation of the moments of a[^^ . One may 
derive or guess some important identities from the expressions for the moments, 
and we also present them. 

In Section |B1 we study various properties of perpetual functionals of the form 

In Section |3 after giving some results on the asymptotic behavior of the 
Laplace transform and the moments of A^^"^ as t tends to oo, we show some limit 
theorems for some Gibbsian measures and moment density measures, following 
Hariya-Yor [H^ . 

In SectionlSl considering another exponential functional of Brownian motion, 
we present a result on a three-dimensional joint distribution. 

Finally, in Appendix ^ we give a proof of an integral representation for the 
density of the Hartman- Watson distribution. 



2. Hartman- Watson distributions 

In this section we give some information about Bessel processes and introduce 
the Hartman- Watson distributions, which were originally obtained and studied 
in from an analytic point of view. 

The starting important fact is the Lamperti relation, which says that there 
exists a Bessel process 

= {Ri^^} with = 1 and index /x satisfying 
exp(B('^)) = i?^^,),, t^O. (2.1) 
For details, see Lamperti 001 and also p. 452 of [57). 
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When /X < 0, we easily see that hnit_^oo ^[^^ — ^oo'' is finite and, from 1)2. Ifl . 
that 

= inf{s; = 0} < oo. 

Dufresne (21] (see also proved that A^^' is distributed as l/27_^, where 
7_yj is a gamma random variable with parameter — /x. We will give a proof in 

Section ini of this article. Baldi et al @j have discussed about in relation to 
some diffusion processes on the hyperbolic plane. 

In view of applications, it may also be worthwhile to note 

where Tc{X) denotes the first hitting time at c G R of a stochastic process X. 
When > 0, replacing Tc by, e.g., last hitting times, we have a similar identity. 

Now denote by pff'^ the probability law of a Bessel process i?^^-* with index 
H on the canonical path space C([0,oo) ^ R) when i?g'^' = x. Then, letting 
TZt = cr{Rs, s ^ t} be the filtration of the coordinate process {Ru}, we have the 
following absolute continuity relationship for the laws of the Bessel processes 
with different indices: for /i ^ and x > 0, 

and, for /j, < and x > 0, 

m..nu...... - ( -Vexpf-^ t7^)-Pi'\.- (2.3) 



^ \ntn{t<ro{R)} J 2 7o (^s)V 

We note that these absolute continuity relationships and an extension of formula 
H2.8|l below also hold for the Wishart processes, which are diffusion processes 
with values in the space of positive definite matrices. For details, see |17| . 

The important relationships (|2.2|l and (|2.3() are shown in the following way. 
Since formula 1)2.2(1 has been presented on p. 450, [Uj and is proven by an easy 
modification of the arguments below, we only give a proof for ((2.3|l in the case 
x = l. 

Let e R, be the probability law of {B^''^} on C([0,cx)) R) and 

{J-t} be the filtration of the coordinate process {^t}. Then, by the Cameron- 
Martin theorem, we have 



= exp(2/iX0M^'"^^ 



We set At{X) = J* cxp{2Xs)ds, a„ = inf{t > 0;At{X) > u} and define 
R = {Ru} by Ru = exp{Xa^)- Then {-/?„, u < to{R)} is a Bessel process with 
index € R. Note that {u < to{R)} = {q:„ < oo} and < W'^'^^Uu < oo) < 1 
if fi < 0. Therefore, by time change, we obtain 

W^''^ - (i?.)2^W^(-^) y^^n{u<MR)} for M > 0, 
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which, combined with (|2.2(1 . imphes 1)2. 

On the other hand, letting p'^'^\t,x,y),t > 0,x,y ^ 0, be the transition 
probabihty density with respect to the Lebesgue measure of the Bessel process 
R^'^\ its exphcit expression is also well known p. 446): ii fj, > —1, 

for cc > and, for a; = 0, 

where is the modified Bessel function. For a proof of (|2.4I) , see Remark 12.11 
below. 

From formula l|2.4|) . assuming i^p'^' = x, we obtain 

E|»p(-A(i!.")=)l.^^^^„p(-^), A>0. (2.5) 

Note that, in fact, in |S7], formula H2.4() has been deduced as a consequence of 

Remark 2.1. As usual, denote by K^^ the other modified Bessel (Macdonald) 
fimction. We set 

Then, by the general theory of the Sturm-Liouville operators, we can show that 
the Green function (with respect to the Lebesgue measure) is given by 

^ e-"*p(^) (t, X, y)dt = 2y ^ F^{V2^x, ^y) 

and, combining this with the following formula for the product of the modified 
Bessel functions (cf. [23|, p. 53) 



I,{x)KM = ^^°°e-«/-(-^+«^)/^%(^)|, < . ^ 



(2.7) 



we obtain p.4() . 

The function i^^ will be used to give an expression for the double Laplace 

transform of the distribution of J^^^ and, when we invert it, formula H2.7|l will 
play a role. For details, see Section^ 

Comparing formulae 1)2.211 and l|2.4|l . we obtain the important relationship 



77 



r ds 



2 Jo (Rs)' 



Rt^y 



(2.8) 
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where Ex^"^ denotes the expectation with respect to Pi°'' . 

The conditional distribution of J^{Rs)^^ds given Rt — y under Pi"\ charac- 
terized by p. 8(1 . is the Hartman- Watson distribution with parameter r = xy/t. 
Of course this is a probabihstic interpretation, and Hartman- Watson have 
introduced a family of probability measures rir{dt), r > 0, on R+, from the 
identity 

For more probabilistic discussions of the Hartman- Watson distributions, see Yor 

m 

In )(j7j an integral representation for the density of rjr (dt) has been obtained. 
Let us define the function 6{r, t), r > 0, t > 0, hy 

0{r,t) = ^^-^e-^/^*^°°e-«V2*e— MC),i,h(C)sin(^)de (2.9) 
Then it holds that Io{r)rir{dt) = 9{r,t)dt, or equivalently 

e-'''*/2^(r,t)dt = /|^|(r), r > 0. (2.10) 

A proof of formula (|2.10l) will be given in Appendix^ We will see in Section^] 
that an expression for the joint probability density of {A[^\ b[^'^) for a fixed t 
may be given by using the function 0(r, t). 

From the explicit expression for 9{r,t), one deduces 

POO 

lim V2^e{r, t) ^ r ^e-''™"''(«) smh{C)d^ 



/•OO 

Jo 



(2.11) 



Remark 2.2. Let {(ft} be a continuous determination of the angular component 
of a complex Brownian motion {Zt} with Zq ~ z ^ 0: 

Zt = |Zt| exp(V^<y5t)- 

Then it holds that 

E[eMV^H^t - ^o)) I \Zt\ =p]= (^i^) . (2.12) 

Moreover it should also be noted that there exists a one-dimensional Brownian 
motion {/?«} independent of {|^s|} such that 



(fit - 1^0= Pat with at = 

Jo 

See also Ito-McKean [37], p.270. 



* ds 



(2.13) 
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Remark 2.3. The function 9{r, t) is, of course, a non-negative function. However, 
since the integral on the right hand side of H2.9|l is an oscillatory one, it is not 
easy to carry out numerical computations, especially when t is small. See 
|36j and |44j for some results and discussions on the numerical computations. 
Remark 2.4. In his fundamental memoir |63| . Stieltjes considered a similar in- 
tegral to that on the right hand side of (|2.9|) . He has shown that 

^ e-^'/^'smh{n£_)sm(^^^d^ = Q (2.14) 

holds for any n G N and for any t > 0. Setting 

/i,,A(dO = ^e-«'/2*(^l -I- Asin(^^))dC, |A| < 1, 
we see from (|2.14|l that the exponential moments of fit,\id£^) are given by 
e"Vt,A(rfe) = e"'*/^ n=l,2,.... 



'R 

and that the moment problem for the log-normal distribution is indeterminate. 
3. Bougerol's identity 

In his discussion of stochastic analysis on hyperbolic spaces, Bougerol has 
obtained an interesting and important identity in law. It plays an important role 
in several domains and in the following discussions. We show it in the simplest 
form. Some extensions and discussions are found in [2] and 0. 

Theorem 3.1. Let {Wt} be a one- dimensional Brownian motion starting from 
0, independent of B. Then, for every fixed t > 0, e:cp{Bs)dWs, sinh(_B4) and 
CLfe identical in law. 

Proof. We follow the proof in |2]. We put 

Xf = exp(Bt) [x + cM-Bs)dWs 
and 

yy = smh{y + Bt). 
Then {Xf } and {Y^} are diffusion processes with the same generator 

1^ ^d^ Id 
2(^ + ")d^ + rd^ 

and, if sinh(?/) = x, {X^} and {Y^} have the same probability law. Hence, we 
obtain the result if we consider the case x — y ^ and if we notice that, for a 
fixed t, Xf is identical in law with e'K.p{Bs)dWs by virtue of the invariance 
of the probability law of Brownian motion under time reversal from a fixed 
time. □ 
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From Bougerol's identity sinh(i?t) WAt, we easily obtain the following 
by considering the corresponding densities and characteristic functions. 

Corollary 3.2. For every a E JH and t > 0, one has 



E 



exp 



1 



exp 



(Argsinh(a))'^ 
2i 



(3.1) 



and 



E 



exp( — -At 



= \— cos(asinh(0)e~«'/2*d^. 

V TTi ./n 



(3.2) 



4. The law of A\ at fixed and at exponential times 

Since the exponential functional appears in several domains and the expressions 
for the density obtained so far are complicated, some authors have tried to ob- 
tain simpler expressions. We give a survey on the results due to Alili-Gruet [3], 
Comtet-Monthus-Yor TS', Dufresne j^H and Yor 69 . For other approaches, see 
Bhattacharya-Thomann-Waymire |HI, Comtet-Monthus 14 , Lyasoff 44|, Mon- 
thus in^l, Schroder jSJ and so on. 

4.1. The first result on an explicit expression for the probability density of 
A\^'^ was given by Yor 69; . 

Theorem 4.1. Fix t > 0. Then, for u > and x £ Ti, it holds that 

P(4'^' e du,Bi'''> e rfo:) =e^--^'*/2^xpf-ii^V(e7u,0^, (4.1) 

\ 2u J u 

where 9 is the function given by (|2.9() . 

Proof. We sketch a proof with fi ~ 0. The general case is easily derived from the 
Cameron-Martin theorem once we have shown this special case. Moreover, since 
a probabilistic proof is given in , which is also found in , we have chosen 
to give a proof based on the general theory of the Sturm-Liouville operators. 

For this purpose, let us consider the Schrodinger operator Hx, A > 0, on R 
with the so-called Liouville potential given by 

Then, by using some properties of the modified Bessel functions, we can show 
that the Green function G(x, y; a^/2) = {Hx+a"^ /2)^^(x,y), a > 0, with respect 
to the Lebesgue measure is given by 

G{x,y;a^/2) = 2/«(Ae^)i^„(Ae^) = 2i^„(Ae% Ae^), x^y, (4.3) 

where the function Fa is given by (|2.6|l . 
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Next we let q\{t, x, y) be the heat kernel of the semigroup exp{—tHx), t > 0. 
Then we have by the Feynman-Kac formula 



qx{t,x,y) = E 



expl-^A^e^-^ 



Bt^y-x 



1 



Jo 



t^y-x) 



~(v-xf/2t 



(4.4) 



'2-Kt 



and 



qx{t,x,y)dt = G{x,y; a'^ /2) . 



Now recall the integral representation (|2.7() for the product of two modified 
Bessel functions. Then, combining the formulae above and setting x = 0, we 
obtain from H2.10|l 



P{At e du\Bt = y) 



cxp[ --U 



la 



2u 

X^{l + e^y)\du 



X^ey\du 
u I u 



2u 



t \dt. 



Using the Fubini theorem, we can easily invert the double Laplace transform 
and obtain 



P{At e du\Bt^y)- 



1 



e-y /2t ^ exp 



1 



=2y 



2u 



e{ey/u,t)- 



du 



which is equivalent to (|4.1() with /.i = 0. 



□ 



Remark 4.1. It may be worth presenting an explicit expression for the heat 
kernel qxit, x,y), which is obtained from 14.4|l : 



Q\it,x,y) 



expl-- 



^ X^ie^"" +e^y)\ fX^e'^+y \d^ 



2C 



.1- 



Another expression will be given by Proposition 15 .41 below. 

Remark 4.2. We find formula H4.3|l in Borodin-Salminen [S], page 264. In that 
book, many explicit formulae related to the exponential functionals are found. 
Moreover, some physicists working with Feynman path integrals have obtained 
similar results in the context of physical problems. We find many formulae and 
related references in Grosche-Steiner |2H|- For example, on page 228 of ^j, we 
find a formula equivalent to (|4.3|l . 

We obtain an expression for the density of Aj''-* from 1)4. However, to 
compare with Dufresne's results mentioned in Subsection 4.3, it is convenient 
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to consider 1/2A['^\ We denote its density by f^^^\a,t). Then, by using the 
integral representation (|2.9|) for the function 6'(r, i), we have 



£ e-«'/2tg-2v^, cosh(c) gi^^i^(^) sin (^^^ 



if ^ > -1. 

When /i = 0, using the elementary formula 

g-VAC^ C Texpf-As-^'j-^^, A,C>0, 

we obtain 



^2(2t)l/2yS7o 



oo 



/(O) (a, i) = ^^—7^ ^ / e-'^t-^'^t"))' cosh(,7)dry 



.-;V2t sinh(20 sin«/t) 

cosh(2e) + cosh(277)^' ^ 



There is also a similar expression for /(^)(a>0 when = 1 (see ^). 

Moreover, by using the formula H4.1|l . we can compute the Laplace transform 
of the function 0(r, t) in r. 

Proposition 4.2. Set adx) — Argcosh(a;) for a; ^ 1. Then it holds that 

for any t > 0. Consequently, one has 

e-^'-d{r,t)dr = -2_a,(^x)a',{x) expl^-^^^ (4.7) 
Proof. Integrating both hand sides of (|4.1() in u yields 



/•oo J 

/ e-"°'^'^(^'e'(r,t) — , 



which is exactly H4.6|l . 

Formula 1)4. 7|l follows from (|4.6I) by differentiating with respect to x. □ 

One finds another proof of Proposition lT!^ in |47j . where the Lipschitz-Hankel 
formula below (cf. Chapter XIII, Watson 6|]) has been used: 

1 

g-rcosh(n) j^^^)^^ = ^ ^ P-^", U > 0, A > -1. 

smh(Mj 
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Several subordinators whose Levy measures are written in terms of the mod- 
ified Bessel functions have been studied in 47 in relation to the exponential 
functional of Brownian bridge. 

Moreover, formula H4.6I) may be proven from the explicit expression 1)2. 9|l for 
the fmiction 9(r,t). In fact, by Fubini's theorem, we obtain 

dr_ _ e^'/^* r e-i'/^' sinh(0 sin«/^) 

r ~ 2V2^ 7r X + cosh(e) 

and, applying residue calculus to the integral on the right hand side, we obtain 
(|4.6|l . By this method, we also obtain for — 1 < x ^ 1 

dr 1 f (Arccos(a;)) 



r°° dr 1 



2t 



4.2. We computed the double Laplace transform of the conditional distribution 

X and obtained an expression for the joint 



of Aj'^-' under the condition 



(m) 



density of {Al''> , Bl'"') in the previous subsection. Alili-Gruet have obtained 
an explicit form of the double Laplace transform of the distribution of A. 



(m) 



Theorem 4.3. For every /i G R, one has 



1 



exp( -^u^A\ 



ip) 



dt = 2 



e^'yFx{u,uey)dy 



for a, u > 0, where A = \J + ^i^ and F\ is the function given by H2.6|l . 

Proof. We use the facts mentioned in Remarks 12.11 and 12.21 At first, by using 
the Cameron-Martin theorem, we rewrite the left hand side as follows: 



exp( ~\u^A. 



E 



t 

exp ( fiBt 
1 



dt 



Wa^ 



exp( ^Bt ~ -u^At 



E 



exp jiB, 



dt 
dt 



u'A, 



IXPt ]dt 



where {Pt} is a Brownian motion independent of B. 

We recall the Lamperti relation and make the change of variables by At = v. 
Then, using the notations in Remark 12.21 and noting that {a^} is the inverse 
process of {At}, we obtain 



E 



exp ^Bt 

= E 



1 



u^At + V^A/3t dt 



Jo 



dv 
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where {Rt} is the radial process of a complex Brownian motion Z starting from 
1 and ifQ = 0. Now, applying formula (|2.12|) . we obtain 



E[{R,Y~^E[e^ 



dt 



dv. 



□ 



Finally, using formulae (|2.4|l and (|2.7I) . we obtain the theorem. 

Monthus has shown another expression for the Laplace transform of the 
probability law of A\'^\ See also |3] for further discussions. 

Moreover, Alili-Gruet have applied the Plancherel identity to the right hand 
side of formula (|3.2|l to obtain 



E 



which has been shown by other methods in and 
Theorem 4.4. One has 

du 1 



P{At e du) = 



VzTTU^ ^/2^^t Jr. 



cosh(^)e^ 



(cosh(C))^/2u-(5+V^7r/2)^/2t 



for any t > and 



lim V2TrtP{A 

t — >oo 



t^zdu)^ -e-^/^'^du, 



u>0. 



To show the theorem, Alili-Gruet 3 have used the Plancherel identity and 
Comtet-Monthus-Yor have used an integral representation for K^/-^^. For 
details, see the original papers. 



4.3. Dufresne [23 has studied in detail the function 



E 



(2A|^Vexp 



2A. 



(4.8) 



which has been shown to be convergent for r e R, Re(s) < 1 and t > 0. He 
has obtained a differential equation in (s,t) which is satisfied by this function, 
explicit expressions in some special cases and a recurrence relation for h^''^{s,t) 
in n. From these results, he has obtained explicit expressions for the density of 

At first we show a computation for the Laplace transform of h^'^'{s,t) in t 
by starting from formula (|4.18|) below, the notation of which we borrow. While 
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this way is converse to Dufresne's original arguments, it is easier to understand. 
For Re(s) ^ and A > ^^/2, one has 



1 



A-/i2/2 
1 



E 



(2A|''V''exp 



2A 



dt 



1 



exp 



\-l?/2T{p) Jo 



where a = (v2A + ^)/2 and /3 = (v2A — /i)/2. After carrying out the integration 
in y and making the change of variables by a; = 1 — u, we can use an integral 
representation for the Gauss hypergeometric function 2F1 (see, e.g., p. 240) 



r(7) 



r(/3)r(7 - /3) Jo 



t/3-i(l _ t)7-/3-i(i _ 7 > /3 > 0, 



and we obtain 



/i'''''(s,A) 



T(a)T(3 + r) , ^ ^ , . 

2il-sr4ia + P+l) + + )■ (4.9) 



After a careful inspection, Dufresne has shown this result under a weaker 
condition on the parameters. 



Theorem 4.5. // v2A > max(— — 2r) and if either (i) Re(s) < 0,r e R, 
or (ii) Rc(s) = 0, r < 1, then formula (|4.9() holds. 

In some special cases, the right hand side of ()4.9|) may be written by using 
elementary functions and we can invert the Laplace transform. 

Corollary 4.6. When r = 1/2 and /i = 0, 1, one has, for Re(s) ^ 0, 

(Argsinh ^/^)^' 



h°^'/^is,t) = 



1 



V2^(l-s) 



exp 



and 



h^^^l^{s,t) - ^exp 
x/2t 



2t 



(Argsinh - 



2t 



(4.10) 



(4.11) 



Remark 4.3. Formula H4.10|l is exactly the same as (|3.1f) when s = —a^. 

We now try to invert the Laplace transform H4.8|l of the probability law 

of (A|''V^- Denote by f'~^'\a,t) the density of (2A|''V^- Then the inversion 
formula is 



/('')(a,i)=e 



2tiJ^ 



e'"'h^'-'\-s,t)ds, 



c 
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where the path of integration C from c — ^/—loo to c + \/—loo with c G R is 
on the right side of the singularities of s, t). Note that the left hand side 

does not depend on r, while the right hand side apparently does. 

When /X = and /i = 1, Dufresne f53 has succeeded in carrying out the 
complex integral and has obtained the following fairly simple expressions for 
/(A')(a,t). 

Proposition 4.7. For t > Q and a > 0, one has 

/(0)(a,i) = ^^= / e-y"^^'e-''^^°''<y^^\oshiy)cos(^]dy (4.12) 
n\'2ta Jo V 2i / 

and 

f(^\a,t) = ^"1,'^' r e-y"/^'e-''^''°'''^y'>^%inh{y)cosHy)sm( ^]dy. 

■nylta Jo V 2i / 

(4.13) 

When /i = 0, we can show the coincidence of the two expressions (14 .Sf) and 
(|4.12|l for /^°)(a,t) by applying residue calculus to the right hand side of 14. 5|) . 
We can also do the same thing when fj, — 1. Combining this coincidence with the 
recurrence formula H4.14|l below also due to Dufresne [53] or H4.15|l . we obtain 
the desired coincidence for general values of /i. 

Theorem 4.8. For all /i, r e R, s with Re(s) < 1 and t > 0, one has 

h^^'^is, i) - (1 - s)^-''/i2'^-^^''(s, t). (4.14) 
Theorem 4.9. For /i e R, a > and t > 0, one has 

7r^/8t-M^t/2 poo 



7r(2a)(''+i)/2Vi 



OO 



^ 7 ^ ) ]Ht^iVasmh{y))dy, 



2 \t 

where Hf^ is the Hermite function with parameter fj,. 

Remark 4.4. We refer the reader to Lebedev 0J| for the Hermite functions. 
We do not need the assertion of Theorem 14.91 in the following and we skip 
explanations. For details of proofs of the theorems above, see the original paper 

m 

Dufresne's recursion relation l|4.14|l has been extended to the level of stochas- 
tic processes and has been written in a more symmetric manner in |50) . The fol- 
lowing theorem plays an important role in studying various diffusion processes 
related to the exponential functional {aI''-*}. For details and a proof, see (SOI 
and [nH. 

Theorem 4.10. Let ^ < v and consider two Brownian motions B^^\ B^^^ with 
different drifts and the corresponding exponential additive functionals {A\.^^}^ 
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{Aj"^}. Then, for every t > and for every non-negative functional F on 
C((0,i] R), it holds that 



F 



A 



< t 



A 



(m) 



F 



A 



+ 2js,s ^ t 



A 



where m = (fi + i^)/2, S — {v — fi)/2 and is a gamma random variable inde- 
pendent of B^'^\ In particular, for every non-negative function ip on R_|_, one 
has 



1 



A[^^J\A\ 



1 



1 



A. 



1 



(^) 



and 



{a-b)^-H"'e''f^''\b,t)db. (4.15) 



4.4. Let Tx be an exponential random variable with parameter A > indepen- 
dent of the Brownian motion B. In and some other articles contained 
in several identities for the distributions of {A^^ , B^j\) ^"^"^ itself have 
been shown. Among them we present two results. 

At first we present an explicit expression for the joint probability density of 



Theorem 4.11. For /i e R and A > 0, set v = \/2X + fi'^ and let p''^\t,x,y) 
be the transition probability density with respect to the Lebesgue measure of a 
Bessel process with index v. Then one has 



P(exp(B^^^) e dy,A)lt' £ du) 



A 



V 



—p^''\u, l,y)dydu. 



Proof. At first we consider the case fj, = 0. Letting /, g : R^ 
negative Borel functions, we have 

<^f\f,g)^E[f{cMBTj)9{ATj] 



(4.16) 



R I be non- 



= XE 



e-^'f{cMBt))g{At)dt 



(4.17) 



Recall the Lamperti relation (|2.1() and, as in Subsection 4.2, put 

du 



as 



(i?i°)) 



which is the inverse process of At. Then, changing the variable in H4.17|l by 
t = a^, we obtain 



<^f\f,g) = \E 



°°e-^«=/(i?(«))5(s) 



(0)^2 
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Moreover, by using H2.8|l and then 12.4|l . we obtain 



which imphes formula H4.16|l when ^ — 0. 

Next we consider the general case. By the Cameron-Martin theorem, we have 

= A / e-^'E[f{exp{Bt))9iAt) eMf^B^ - fih/2)]dt 
Jo 

/>oo 

= A / e-''"'/^E[f{eMBt))9iAt)eMf^Bt)]dt 

JQ 

2A 

= —E[f{exp{Bs,^)9{As^)eMf^BsJ], 

where Sn (= T^i ji with our previous notation) is an exponential random variable 
with parameter j1 independent of B. 

We apply H4.16|l in the case /i = 0, which we have just shown, to obtain 



C30 /'C30 



(/, 5) = A / / g{u)j{y)y-'^''^^-^^v^^^ {n, 1, y)dudy, 



'0 "'0 

which implies formula H4.16|l in full generality. □ 
The following result was used in the previous subsection. 



Theorem 4.12. For /i e R and X > 0, set v = a/2A + /x^, a = + n)/2, 
b = (ly — fi)/2 and let Zi^a and be o, beta variable with parameters (l,a) and 
a gamma variable with parameter b, respectively; 



P{Zi^a e dt) = a(l - t)°"^dt, < t < 1, 
P(76 e dt) = j^t^'-^e-'dt, t > 0. 

r(6) 



Then one has the identity in law 



where, on the right hand side, Zi,a and «re assumed to be independent. 
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Proof. Recall the Lamperti relation H2.1|l and use the same notation. Then we 
obtain 



P{Ai^^ ^ m) = A / e-^'PiA^"^ ^ u)dt 

poo 

= A / e-^*P(a(r) ^ t)dt = E[cxp{-Xai^^)], 



where 







al^) - / , , dv. 

Jo {Ri^^y 



Moreover, by the absolute continuity relationship H2.2() or (|2.3() for the laws of 
different dimensional Bessel processes, we obtain 

where P^,"' = r'^q^ = 1. 

Then, recalling formula H2.5I) and using the elementary formula 

/•oo 

/ e-^^'x''-^d\ = T{b)C^'', 
Jo 

we obtain after a simple change of variables 

1 r°° 



1 



l/2u j-l 

v^'^e'^dv / ail-ty-^dt, 

J2uv 



m 

which shows H4.18|l . □ 
5. Moments 

In this section we present some results on the moments of A[^^ and their con- 
sequences. We show the results on the positive moments in detail and we only 
consider negative moments in some special cases. We note that the expecta- 
tion £'[log(yl|''^)] has been also considered in the context of physics in Comtet- 
Monthus-Yor (15) . 

5.1. In this subsection we present several results about positive moments of At. 
At first we present an easy consequence of the Bougerol identity sinh(i?t) = 

Proposition 5.1. For any t > 0, one has 

EU,n = /^(sinh(.))-^e-^V.*,,, . = 1,2,... 
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From the proposition, we obtain 

E\(Atr] = — , , e^"'*(l + o(l)) as t ^ oo. 

' ^ r(n + 1/2)23"-! V V " 

To compute £'[(yl|'''')"] for a general value of /i, we note 

E[{A[^^r] = e-^'*/2^[exp(/ii?,)(A)"], 

which is easily shown by the Cameron-Martin theorem. For the expectation on 
the right hand side, we have the following (cf. |35j.|65j .16911. 

Theorem 5.2. (i) Set (j){z) — Then one has 

/>oo n, ^ 

I e--i^[exp(Mi^O(^0"]^^-»!n ,,^(^^,^.) (5.1) 

for any n G N and /i ^ i/Re(A) > + 2n). 
(ii) For any i > 0, n G N, /i ^ 0, it holds that 

n 

E[exp{^iBt){Atr] =Y,C{n,j;fi)eMt<l>{l^ + m, (5-2) 

where C(0, 0; yit) = 1 and for j — 1, n 

C(n,j;/.) = 2-"(-l)"-^Y") n {^, + J + k)-\ 



0<k<7 



Proof, (i) By the independent increments property of Brownian motion, we 
deduce 



■ nlE 



exp(^Bi) ■■■ exp(2(Bsi H h BsJ)dsi ■ ■ ■ dsr 

• • • / E[eM^l{Bt - BsJ + + 2)(B,„ - B,„_J 

+ • • • + (^t + 2n)i?si)]dsi • • • dsn 
■ ■ ■ exp{(j){fi){t - s„) + (f>{^i + 2)(s„ - s„_i) 

+ ■ ■ ■ + 0(m + 2n)si)(isi • • • dsn- 
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Hence, setting sq = and s„+i = t, we obtain 

/>oo 

Jo 



i=o 

(ii) By the additive decomposition formula, we have 

TT i = y ^ 



where 



Since 



t=a ^ - -^(m + 2j) A - <^(/x + 2j) ^ 



^r- n 



0</c<n 



1 



we can invert the Laplace transform (|5.1|l and, after some elementary compu- 
tations, we obtain (|5.2|) . □ 

The same arguments show that formulae (|5.1|l and H5.2() hold when is any 
complex number. In particular, for /.j = \/—la, we can show 

n 

E[exp{v^aBt){Atr] = ^ C(n, j; \/^Q)e'^(^"+2^)*, a G R. (5.3) 

By Fourier analysis, we obtain the following result on the conditional mo- 
ments from (|5.3|) . For details, see j35) . 

Proposition 5.3. For every t > and n E N, one has 



1 



2TTt , , 

(5-4) 

5e-^ /2*(cosh(fe) -cosh(|a;|))"d6. 



Recall that the Bessel function Jo is defined by 

■^o(^)-Z. ■ 
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Then, at least heuristically, with the help of this series expansion and inter- 
changing the order of the integration and the infinite sum on the right hand 
side of (|5.5|l . the next formula holds: 



In fact, this identity has been shown in |H] and |SSI by different methods. 

Comparing IjS.Sfl with H4.4|l . we obtain the following expression for the heat 
kernel q\{t,x,y) of the semi-group generated by the Schrodinger operator H\, 
A e R, with the Liouville potential given by H4.2|l . 

Proposition 5.4. For t > 0,x,y E H, it holds that 



Remark 5.1. From formula 1)5. 6|l . we can show the well known explicit expression 
for the heat kernel of the semigroup generated by the Laplace-Beltrami operator 
on the hyperbolic plane. For details, see \W\ and jSS]- The Brownian motions 
on hyperbolic spaces will be discussed in the second part of the survey |^ . See 
also and jJTj about the close relationship between the Brownian motions 
on hyperbolic spaces and our exponential functionals. 

Remark 5.2. Denoting by C„ the right hand side of formula H5.4|l . it is easy to 
show 



This means that Carleman's sufficient condition for the unique solvability of the 
Stieltjes moment problem does not hold for the (conditional) distribution of 
Recently Horfelt |34) announced that he has proven the indeterminacy for the 
moment problem for At- 

Dufresne has shown that the Stieltjes moment problem for the reciprocal 
(^1'^'')"^ is determinate. By using this fact and formula H4.1|l . he has also shown 
another representation for the density of ^4^''^ after giving another expression 
for the moments. See also for a similar result for exponential functionals of 
Levy processes. 

We go back to formula (|5.4(l . Let us consider an exponential random variable 
e whose density is e~^, x ^ 0. Then we have ii^[e"] — nl and 




(5.5) 




(5.6) 



OO 




b > \x\. 



Hence, assuming that e and B are independent, we obtain from (|5.4|) 
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Although we do not know the unique solvabiUty of the moment problem for the 
conditional distribution of At, we can show the following result. It says that, if 
we multiply At by an independent exponential random variable, the probability 
law of At is modified into a rather simple one. 

Theorem 5.5. For any t > 0, the conditional distribution of eAt conditionally 
on Bt — X coincides with that of e^(cosh(V2et) — cosh(a;)) given e > x^ /2t. 

For a proof, see PHI or gS] . 

In these articles one also finds the following results, which are shown by using 
the expression (|4.1|) for the joint density of {At, Bt). 

Theorem 5.6. (i) For any t > and A > 0, one has 



E 



exp 



Bt 



exp 



2t 



(5.7) 



where ifx{^) = Argcosh(Ae~^ + cosh(x)). 

(ii) Letting e be an exponential random variable independent of B and {Lu,u ^ 
0} be the local time of B at 0, one has the identity in law for fixed t 



{eAt,Bt 



(law) 



(e^'(cosh(|Bt| + Lt) - cosh(Bt)), Bt). 



Proof. We only give a proof of formula l|5.7|l . which will be used in the next 
subsection. Let / be a non-negative Borel function on R. Then, from formula 
(|4.1|l . we deduce 



f{^)E 



R 

= E 



exp 



Bt=x 



1 



f{Bt)cxp[- — 



f{x)dx I exp 
R Jo 



A + cosh(a;) 



0^e.yu,t) — 
u 



f{x)dx / exp(-v3:r(A)w)6'(w, t) 



dw 



Then, with the help of formula (|4.6() . we obtain 



I{^)E 



R 



expl-^ 



Bt^x 



1 



f{x)^= 

R^ '^t 



exp 



^.(A)2 



2t 



e--"'^'dx 
dx. 



□ 



5.2. As is mentioned in Remark 15.21 it is important to have formulae for the 
negative moments for A\.^\ Moreover, considering the negative (conditional) 
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moments of A[^\ we obtain several interesting identities which are related to 
some special functions. However, we restrict ourselves to some special cases 
which can be obtained from the results mentioned in the previous sections. For 
more results, see [B], HO], [213, [221 ^^'^ the references cited therein. 

We again start from the Bougerol identity or rather from formula H3.1|l . 



Proposition 5.7. For any t > and a £ R, one has 



wl/21 ^ ^-1/2 



id E 



cxp{aBu)du 



^1/2-1 



= 1. 



The first identity is a direct consequence of formula (|3.1() . and, from this, we 
obtain the second striking one by the scaling property of Brownian motion. 
The following is also a consequence of the Bougerol identity. 



Proposition 5.8. For any p with Q < p < 1/2, one has 

r(a) ^-1/2, 



r'/^(l + o(l)) as t 



oo. 



(5.8) 



In fact, Hariya-Yor j321 have shown that formula H5.8|l holds for any p > 0. In 
1221 one finds a complete description on the asymptotic behavior of E[{A['^^)p], 
p G R, as t — > oo and on the limit of the probability measure 

{E[iA[^^r]r\A[^^rdP. 

We will present a part of the results of ^2 in Section [T] 

To find an analogous result to that of Proposition 15.71 for the conditional 
distribution of At given Bt, we recall formula H5.7|l . By differentiating both 
hand sides of (|5.7(l in A, we can compute E[{At)~"-\Bt = x], n = 1,2,... In 
particular, we obtain the following. 

Proposition 5.9. For any t > and x £ R, one has 



E[iAt)-'\Bt ^ x] = 



2d E[{At)-'\Bt = 0] = t- 



t sinh(a;) 

By the scaling property of Brownian motion, we also obtain 
f-i 



E 



exp(Q:i3„)c!u 



Bi = 



= 1 



for any a G R. 



Chaumont-Hobson-Yor 13 have shown that this striking identity is a conse- 
quence of the so-called cyclic exchangeability property of the Brownian bridge 
and they have obtained several similar identities for the stochastic processes 
which enjoy this property. 
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6. Perpetual functionals 

Throughout this section we assume > 0. For a non-negative Borel measurable 
function /, we set 



Jo 



and caU it a perpetual functional. 

At first we show a necessary and sufficient condition in order that A^"^ (/) is 
finite almost surely. We follow Salminen-Yor [35]. See also Engelbert-Senf pi| 
for another approach. 

Theorem 6.1. Let ^ be positive and f be a non-negative and locally integrable 
Borel function on R. Then, in order for (/) to be finite almost surely, it is 
necessary and sufficient that J f{y)dy < oo. 

Before proceeding to a proof, we recall a lemma due to Jeulin j88| . 

Lemma 6.1. Let {Zy,y ^ 0} be a stochastic process such that the probability 
law of Zy does not depend on y and is absolutely continuous with respect to 
the Lebesgue measure. Assume also that E[\Zy\] < oo. Then, for any a-finite 
non-negative measure v on R+, \Zy\v{dy) is finite almost surely if and only 
if i^(R+) < oo. 

Proof of Theorem \6.1\ Letting L^f^^{t,y) be the local time of B^f^^ at y up to 
time t, we have 



f f{Bi^^))ds^ f f{y)L^^\t,y)dy 
Jo Jr 



and, letting t tend to oo and using the monotone convergence theorem, we obtain 

A^£\f)= I f{y)L^-^K^,y)dy. 
Jr 

The distribution of L^^-* (oo, y) is known (cf. 59 ) and is given by 

PiL^^\^,y)^v)^hZ:, 



y<0. 



The probability law of {i'^^^(oo, y), y ^ 0} does not depend on y. Hence, 
letting / be the non-negative function given by f(x) — f{x)l-R^{x), we obtain 
from Lemma l6. II 

p(^J^ f{Bi^''>)ds < oo^ = 1 if and only if J f{x)dx < oo. 
The rest of the proof is easy. □ 

In [SHI, sufficient conditions on / for A^\f) to have moments of any order 
and to have exponential moments are also given. 
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We next consider the important and typical case, f{x) = e We set 



exp(2B(-''')d5 



(law) 



exp 



The following Theorem 16. 21 is due to Dufresne |21|. 

Bertoin-Yor 0| (see also 0) have considered Levy processes which drift to oo 
as i — > oo and have shown that the distributions of the corresponding perpetual 
functionals are determined by their negative moments. Moreover, they have 
given another proof of the following theorem. See also \57\, p. 452, and |7U| . 

Theorem 6.2. Let fi > 0. Then, A^^^-* is distributed as {2jfj_)^^, where 
denotes a gamma random variable with parameter fi. 

Proof. Differently from the proofs given in the above cited references and at the 
end of this section, we give an analytic proof. 



Let {ei'^\t)} be the diffusion process given by ex^'{t) = xe'xp{B]:^') and T; 
be its first hitting time of z. For a > 0, we set 



Vz{x; a) ^ E 

Then Vz solves the equation 
1 



cxp(-y / ie(d'\t))-'dt 



v{z) = 1, 

Moreover, we have 

lim Vz{x; a) 



limwfa;) = 0. 



exp 
exp 



a 

2a;2 °° 



-2B\ 



dt 



(6.1) 



Since Vz{x] a) depends only on x/a, we denote it by Vz{x/a). 

On the other hand, making the change of variables ^ = a/x and Vz{x/a) 
^6 obtain after an elementary computation 



This is the differential equation for the modified Bessel functions and of 
order /i. Since 0(^) — > as ^ — s- 00 by the boundary condition for u^, we easily 
obtain 



a J Kf^{a/z) 
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xj Kf,{a/z) 



and, therefore, 

Vz{x/a) 
Using the formula 

K^{l^)=2^-^T{^i)1^-^'{l + o{\)) as 77^0 
and setting x = 1, we get from (I6.1() 



E 



exp 



1 



r(M) V2 



(6.2) 



Finally, with the help of the integral representation for (cf. ^41,, p. 119) 



1 /a 



^-(«) = 2 2 



□ 



and the injectivity of the Laplace transform, we obtain the assertion. 

Theorem 16.21 has been extended in several ways. For example, in j58j . the 
following has been shown. 



Proposition 6.3. Setting 

one has 
E 



f 

Jo 



e^^^^^'l , ^ ds 



exp(^-i(a2^(^^-+)+/32A(^A<-)) 



(6.3) 



r(/i)2A'-i (3K^{a)I^-i{P) + aK^_iia)I^{P) 

for any a, /3 > 0. 

We recover formula (|6.2() from 16.3|l by taking a = P and using the formula 
i^^(z)Vi(z)+K^_i(z)/^(z) = i. 

Moreover, by the Lamperti relation, we have 

A(^^)=infK4-^)=0}, 

where is a Bessel process with index —fi<0 starting from 1. Hence, 

by standard diffusion theory, we can verify 



E 



exp 



1 



2 /a 
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and we have another proof for Theorem 16.21 

In Salminen-Yor 1201 , the distributions of several perpetual functionals have 
been identified with those of the hitting and occupation times of Brownian 
motion and Bessel processes. See also ^Hl for more results. 

When /i ^ 0, v4|^'' oo as t ^ cx3. We now show the following limit theorem 
forlog(4'^)). 

Proposition 6.4. Let Z be a standard normal random variable. Then, as t — > 
oo : 

(i) = 0, t^^/'^ log{A^f'^) converges in law to 2\Z\; 

(ii) > 0, i"i/2(iog(^(M)-) _ 2^t) converges in law to 2Z. 

Proof, (i) By the scaling property of Brownian motion, we have 

,(0) (la*') / 2s/tB, 



and, by the Laplace method. 



■^\og( [ e^^^^ds^ ^ 2 max Bs. 



Since maxo<5<i Bg — \Z\, we obtain the result. Note that another proof may 
easily be obtained with the help of Bougerol's identity. 

(ii) By the invariance of the probability law of Brownian motion by time reversal 
from a fixed time, we have 

hence 



Since t~^/^Bt is identical in law with Z and the second term on the right hand 
side tends to as < ^ oo almost surely, we obtain the result. □ 

By using the Markov property of Brownian motion, we can show the following 
process level extension of Proposition 

Proposition 6.5. In the sense of convergence of finite dimensional distribu- 
tions, 

(i) if fj. = 0, {X^^/'^log{A^^^),t > 0} converges to {2 maxQ< Bs,t > 0} as 
A oo ; 

^ n r \-i/2n„„/' /^^^ 



(ii) > 0, {A-i/2(-iog(^(M)) _ 2^Xt),t> 0} converges to {2But > 0}. 



We end this section by presenting some results on the asymptotics of the 
joint distribution of Bt and At as t — > oo. 
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In Section El we showed 



hm V2^6'(r, t) = Ko{r), r > 0. 

t — >oo 

From this result, we can show the following. 

Proposition 6.6. For any non-negative functions f onR and g on R+, satis- 
fying 

\x\f{x)dx < oo, / f{x)dx < oo, / g{y)e~'^^^^ — < oo, 







\imV2^E[fiBt)giAt)]^ [ H f{x)g{y)K^{e- ly)e^J- 



one has 

1 + e^^ ^ dxdy 

This result has been extended by Kotani and, recently, by Hariya [3J to 
general Brownian functionals of the form V{Bs)ds for a non-negative func- 
tion V, instead of the exponential function. They have shown, under suitable 
assumptions, that 

rt 

lim t^^^E 



/(St).9^ ViBs)ds 



exists and may be identified. 

Kotani used Krein's theory of strings and Hariya obtained the result by re- 
ducing the computations to those for the three-dimensional Bessel process. For 
details, see the original papers. 

Remark 6.1. As is mentioned above, the exponential perpetual functionals of 
Levy processes have been studied in [5] and [7]- 

On the other hand, for two independent Levy processes {Yf} and {Zt}, which 
are given by the sums of constant drifts, Brownian motions and compound 
Poisson processes, the perpetual functional exp{Yt)dZt has been studied in 
Gjessing-Paulsen [22], Nilsen-Paulsen Paulsen [2^1 and so on. For details, 
see the original papers. 

7. Some limiting distributions 

In this section, following Hariya- Yor [^, we present the results on the asymp- 
totic behavior of 

i^^"' = £;[exp(-a4''^)] and A,^''''") = 

as t — > oo and on the weak limits of the probability measures 

p^''-" V exp(-a4''V-P and (A*'''" 

We present the results without proofs. For details and proofs, see the original 
paper. 
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Theorem 7.1. Fix a > 0. Then we have the following: 

(i) // > 0, it holds that 

(ii) If fi = 0, it holds that 

lim tV2^(o,a) ^ [lKo{s/2a). 

(iii) If ij. < 0^ it holds that 

t-oo * r(|/i|) (v^)/* 

For a|^'™\ we consider the expectation of a more general functional. For 
a ^ 0, ^ > 0, we set 

To present the result, we need to partition the (/i, m)-plane into the following 
six subsets: 

i?i = {2m > /i, ^ > 0}, Li = {2m = fj,, ^ > 0}, 
i?2 = {^Tt < /i, 2m < /i}, ^2 = {m = /i, /i < 0}, 
Rs = {m > n, 0}, Lg = {m > 0, = 0}. 

Here are some corresponding limiting results. 
Theorem 7.2. (i) If{fi,m) € -Ri, one has 

limt3/2e/'^V2A('''™)(a,0 

2(,.-5)/2 / X / \ .oo ^m-l-^/2g-a:/2 



(ii) // (/i, m) e -Li , one /las 



hm ti/2e^ */2a(^'")(o,0 = ^-^r(TO)r 



(iii) If {ii,m) e i?2, one /las 



lim e2™(''-™)*A^'''"V>0 = 



2"T(^-m) 
T{ii-2m) ^ 



H. Matsumoto, M. Yor/ Exponential functionals of BM, I 



339 



(iv) If {iJL,m) e L2, one has 



lim t-iA^'")(a,e) = ^^r 



(v) // {p,m) e i?3, one has 



r>fi /-oo m-fj,-l -x/2 

lim A ^'"^(a,0 - 7^ / ^7 % dx. 



(vi) If {^,m) e L3, one has 



Next, letting yV'(^) be the probability law of B*^^^ on the canonical path 
space W — C([0, 00) R), we consider the Gibbsian-like measure pj^'"-* and 
the moment density measure p^^^'™^ on W defined by 

df[^-'^\x) = {d[''''^^)-^ eyi^{-aAt{X))dW'^^'\X) 

and 

for /i e R, a > and m e R, where At{X) — J* exp{2Xs)ds. 

One of the aims in j32j is to show that the probability measures p^'''"-' and 
-p^(P:™) converge weakly on as t — > 00 and to characterize the limiting mea- 
sures. 

To present the results, we need to recall the generalized inverse Gaussian 
distribution, which we denote by GIG(/i;a, 6), and its relation to Brownian 
motion with drift. 

For a, 6 > and /i G R, the generalized inverse Gaussian distribution is the 
probability distribution on R_|_ whose density is given by 



-\ b^x-\ ]]dx, X >0 



In particular, GIG(— 1/2; a, 6) is called the inverse Gaussian distribution. 

Let fj, > and denote by ri^y and ii'^y the first hitting time at y and the last 
exit time from y of {x + B^^^^t)}, respectively. Then, it is well known that the 

distribution of ri^y is GIG (—1/2; y—x, fi) and that of ii'^Xy is GIG (1/2; y~x, 11). 
When X = y, iy^y is distributed as 2/^~^7i/2 for a gamma random variable 71/2 
with parameter 1/2. By the strong Markov property of Brownian motion, we 
have for x < y 
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which yields the special case of the following identity in law ()7.1|1 (with 1/ — 1/2, 
a = y — X, b = 11). 

Generally, we let ij^^""^ be GIG(±i^; a, b) random variables and ji, be a gamma 
variable with parameter u. Then we have 



(,) (law) . 



for all > 0, where we assume that the two random variables on the right 
hand side are independent. We can easily prove this identity by computing the 
Laplace transforms. For more details of this identity, see 021 ^md |62| . 

For our purpose, we only need to discuss a particular case and let J(/i, a) be 
a GIG(/i; a, 1) random variable. 

Moreover, letting T^, z > 0, be the transformation on the path space W given 

by 

T,{X){t)^Xt-log{l + zAt{Xj) 

and assuming that I{^,a) is independent of B^^^\ we denote by W'j|'^ the 
probability law of the stochastic process {r/(-^ „)(_B'^^))(t), i ^ 0}: 

Theorem 7.3. (i) // /i ^ 0, then p|^'"^ converges weakly to wj"p ^2") ""^ ^ 
as t oo. 

(ii) // /i < 0, then pj'*'"-' converges weakly to ^2")' '■^ ^Qual to 

Note that, if /i ^ 0, the limiting measure does not depend on /i. Some heuristic 
discussions and analytical explanations are found in the original paper j32j. 
For the moment density measures "p^f^'"^^ ^ we have the following results. 

Theorem 7.4. Let 7^ be a gamma random variable with parameter v > 
independent of B^f^\ Then, the following holds: 

(i) If {^,m) G i?i U L3, 'pj:'^'"^^ converges weakly to 14^27^ /2 on W as t —> 00; 

(ii) // (At, m) e Li U i?2 U L2, v[^'"'^ converges weakly to >v(M-2m). 

(iii) If {^,m) G i?3, 7?^^^'™-* converges weakly to W2~^'f'_^. 

Remark 7.1. Both the generalized inverse Gaussian distributions and the path 
transformations will also play important roles when we study some diffusion 
processes associated with the exponential functionals. For details, see [THl ISU] . 
This will be discussed in the second part of our survey |51| . 



8. A 3£) distribution 



In the different studies of the so-called stochastic volatility models in mathe- 
matical finance or in those of generalized Laplacians on the hyperbolic plane 
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and usual Laplacians on complex hyperbolic spaces, or again in those of some 
diffusion processes in random environments (see, e.g., ^Hli lUli |32]> EHI)) 
one often needs to consider, together with A['^\ another exponential functional 



exp 



and to study the joint distribution of {a[^\ 

In this section, following 1,, we show that our method to prove Theorem 14. II 
may also be applied to this study and that we can give an explicit expression for 
the Laplace transform of the conditional distribution of A\.^'^ given {af"^ , b[^^ ) . 

To present the result, we denote by i/jj'^-' {v, x) the joint density of (a^^'' , i?^^'' ) 
given by 



1 

16^ 



'v exp 



which is easily shown from (|4.1() by the scaling property of Brownian motion. 
We also set 

2Aexp(x/2) 
(j) = (f>{v, x; A = . , ■ 
smh(Aw/2) 

Theorem 8.1. For any t > 0, X > 0,v > and x e R, one has 



E 



Am) 



exp(-A(l + e^) coth(Aw/2))6i(0, t/4). 



4'[^\v,x) 



A 



4sinh(Aw/2 



Proof. We only give a proof when /i = 0. Once this is done, the result for 
a general value of fj, is easily obtained by the Cameron-Martin theorem. We 



denote A^^^ and a^^' simply by At and a^, respectively. 

For A > and fc e R, we consider the following Schrodinger operator H\ 
with the Morse potential V\^k on R: 



2da;2 



Letting q\^k{t, x,y) be the heat kernel of the semigroup generated by H\ k, we 
have by the Fcynman-Kac formula 



q\,k{t,x,y) 



= E 



exp 



-^A^e^-At + ^Afce^at 



Bt^y-x 



^-\y-x\'/2t 



The Green function G\_k for H\ k is obtained by using the general theory of 
the Sturm-Liouville operators. In fact, letting Wk,a and M^^a be the Whittaker 
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functions (cf. ^2], El) 7 can show 

/>oo 

Gx,k{x, y; o? ji) = / e-''"'/\x,k{t, x, y)dt 







Ai (1 + 2a) 

for y ^ X and a > 0, where, if A: > 0, we assume a > k — 1/2. See also jHj p. 196. 

Note the integral representation for the product of the Whittaker functions 
(see Hg, p.86, and [211, p.729): 



Then we obtain 

Jo sinh(u) \ smh(u) / 

by a simple change of variables, e" ~ coth(z/2). 

Recalling the characterization 12.1U|I for the function d{r,t), we obtain 

roo 2ku 

Jo 2smh(u) 

by the injectivity of Laplace transform, where = 2Ae'^^+^'/^/ sinh(u). The rest 
of proof is easy and is omitted. □ 

Another probabilistic proof of the theorem is given in 

Appendix A: Density of the Hartman- Watson distributions 

In this section, following [HZI , we show that the density of the Hartman- Watson 
distribution is given by {Io{r))~^6(r,t), where the function 0{r,t) is given by 

Theorem A.l. Let 9{r, ■),r > 0, be the function on (0, oo) given by (12. 9|) . Then 
we have 



e-''^t/^e{r, t)dt = I^{r), > 0, r > 0. 

Proof. We start from the following well known integral representation for the 
modified Bessel function Iu{z) (see, e.g., jHl], p. 176): 

h{z)^ \= I e^'=°''^("')-'^'"(i«;, Re(z) > 0,Re(z^) > -1/2, 

27rv — 1 Jr. 
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where the path C of integration consists of three sides of the infinite rectangle 
with vertices at oo — tt\/—1, —Tr\/—1, tt^/—1 and oo + ti^/—\. Using the elemen- 
tary formula 



V27rt3 



we have 



e - ■''^e{r,t)dt, 



27^^/^ 



c 



We write down the complex integral as the sum of the integrals along the half 
lines and the line segment to obtain 



-1 



27rV-l Jo 
1 



1 

2W^ 



The second term on the right hand side is zero because the integrand is an odd 
function. Hence, we obtain after a simple computation 



TT 

rt 



-rcosh(f) -(flit, 
e ^^'e ^ ' TT COS 



,^V2t 



-T cosh(^) , 



sinh(Oe-« /^tgjj^i :ii j^^ 



— ^ sin 



d^, 



t 



where we have used integration by parts formula for the second equality. 

Now, by the injectivity of Laplace transform, it is easy to show 6{r, t) = 
e{r,t). □ 
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